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Abstract 

A procedure is presented for efficient generation of 
high-quality unstructured grids suitable for CFD 
simulation of high Reynolds number viscous flow 
fields. Layers of anisotropic elements are generated by 
advancing along prescribed normals from solid 
boundaries. The points are generated such that either 
pentahedral or tetrahedral elements with an implied 
connectivity can be be directly recovered. As points are 
generated they are temporarily attached to a volume 
triangulation of the boundary points. This triangulation 
allows efficient local search algorithms to be used when 
checking merging layers. The existing 
advancing-front/local-reconnection procedure is used 
to generate isotropic elements outside of the anisotropic 
region. Results are presented for a variety of 
applications. The results demonstrate that high-quality 
anisotropic unstructured grids can be efficiently and 
consistently generated for complex configurations. 

Introduction 

Unstructured grid technology is a promising approach 
offering geometric flexibility for handling of both 
complex geometry and physics. As such, it can provide 
a powerful capability for accurately and efficiently 
computing complex flow fields about realistic 
aerospace configurations. Several unstructured grid 
generation and flow solver procedures have been 
developed and successfully demonstrated for inviscid 
flow about complex configurations. For isotropic 
elements, existing procedures are robust and capable of 
generating high-quality grids efficiently. For 
anisotropic elements in viscous flow applications, 
further improvements in efficiency, robustness, and 
quality of unstructured grid generation procedures are 
needed. In this paper, a grid generation procedure is 
presented which offers the potential for overall 
improved performance and quality. 

The most common approach used to generate 
anisotropic unstructured grids is to use a layered 
approach and generate points along normals from solid 
boundaries. Unstructured grid generation for viscous 
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applications have been developed using a modified 
advancing-front method by Hassan, el al*, a 
semi-structured approach by Lohner^, 
advancing-normal by Marcum-^, and advancing-layers 
by Pirzadeh"^. Hybrid methods for prismatic/tetrahedral 
grid generation have been developed by Kallinderis*^ 
and Sharov and Nakahasi^. While all of these methods 
differ in how points and elements are generated, they all 
produce very structured and aligned elements adjacent 
to solid boundaries and use isotropic tetrahedral 
elements outside of the anisotropic or boundary-layer 
region. Use of prismatic elements within the anisotropic 
region can reduce subsequent memory and CPU 
requirements for the flow solver without any loss of 
accuracy. 

The goal of the present work is to modify the existing 
advancing-normal and advancing-front 

local-reconnection method^ for efficient generation of 
high-quality, mixed element type unstructured grids for 
viscous flow applications. While this method has many 
advantages, the local-reconnection process is unable to 
remove trapped sliver elements. Such elements can be 
formed between prismatic groups of elements as shown 
in Fig. 1. Local-reconnection or connectivity 
optimization can not remove these trapped slivers as 
they represent a local minimum state that cannot be 
removed without reconnecting a potentially very large 
number of elements. An alternative is to use the same 
point placement strategy and discard the connectivity in 
favor of a hybrid approach. With a hybrid approach the 
element connectivity is directly implied. Also, the 
elements can be recovered as either all tetrahedra or a 
mixture of five and six node pentahedra and tetrahedra. 
This combined approach retains most of the generality 
of the original and is very efficient. 

Unstructured Grid Generation Procedure 

The approach used in the present work uses the 
advancing-normal point placement algorithm^ to 
generate points within the anisotropic region. The 
advancing-front/ local -reconnection (AFLR) 
procedure^-® is used to generate tetrahedral elements in 
the isotropic region. The basic steps in the overall 
procedure are listed below. 

1 ) Generate a boundary surface grid. 

2) Generate a volume triangulation of the boundary 
points. No boundary recovery is required. 
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3) Create new points using advancing-normal point 
placement. 

4) Attach new points to the existing triangulation for 
searching and checking. 

5) Create a new boundary surface grid using the 
inflated surface from advancing-normal point 
placement, 

6) Use AFLR to generate an isotropic tetrahedral 
element grid for the remaining regions. 

7) Merge anisotropic and isotropic regions. Element 
connectivity within the anisotropic region is directly 
determined from the point ordering. 

Advanc ing-Norm al Anisotropic Grid Gen e ration 

With advancing-normal type point placement for 
high-aspect-ratio elements, the standard AFLR 
procedure^ does produce sliver elements of the type 
shown in Fig. 1. These elements are generated only in 
regions of high-aspect-ratio elements with a very 
structured alignment. Elimination of these elements 
with local-reconnection is not feasible. There may be 
no nearby optimization path which produces a better 
connectivity. The problem is inherently due to the very 
structured nature of the grid in these regions. Only a 
limited set of possible triangulaiions, that do not contain 
sliver elements, exists for a set of tetrahedra aligned in 
prismatic groups. A modified process is proposed here 
which eliminates the sliver problem and retains the 
generality and efficiency of the original procedure. In 
the present approach, local-reconnection is not used to 
determine the connectivity in these regions, Instead, the 
connectivity is directly determined by the order in 
which points were generated. This produces a very 
structured connectivity and allows the elements created 
to either be all tetrahedra or of mixed type. 

I 

Fig. 1. Trapped sliver element between prismatic 
groups of tetrahedral elements. 

The basic steps in the modified procedure are listed 
below. 

1) Determine a normal vector at each active 
boundary-layer point. Initially the normals arc 


based solely on the original boundary surface 
geometry. As the generation advances the normals 
are generated using the geometry of the outer layer 
of the boundary-layer grid. 

2) Smooth the normal vectors with a weighting 
dependent upon the distance from the boundary. 
Initially the normals arc unsmoothed. At the 
estimated end of the boundary-layer region full 
smoothing is applied. The end of the 
boundary-layer region is based on a estimate of 
where the element aspect ratio will be near isotropic. 

3) Generate new points one layer at a time. New points 
are created along the normal vector with the normal 
spacing determined using geometric growth from 
the boundary surface. Generation of new points 
along a normal is shown in Fig. 2. 

4) Check distance between new points and surrounding 
element quality. The volume triangulation is used to 
efficiently check nearby points. As boundary-layers 
merge new points may be too close and 
advancement should terminate locally. A new point 
is rejected if the distance between it and any nearby 
new (or existing) point is less than a preset fraction 
of the local element length scale. Boundary-layer 
advancement is terminated locally if a new point is 
rejected. 

5) New points arc also rejected if any of the 
surrounding elements that they may produce fail a 
quality check (maximum angle < 160 deg.). 
Boundary-layer advancement is terminated locally 
if a new point is rejected for quality. 

6) Active points can become isolated as 

boundary-layer advancement is locally terminated. 
This can be prevented by rejecting a new point and 
terminating local advancement if more than some 
fraction of its neighbors have been terminated. 

7) Check element aspect-ratio. As the grid advances 
and the normal spacing increases the element 
aspect-ratio will eventually be isotropic. 
Boundary-layer advancement is terminated locally 
when the aspect-ratio on the next layer would be 
greater than unity. 

8) Attach accepted new points to the volume 
Iriangulation. New points are connected and 
attached to the existing element that contains them. 

9) Generate a new boundary surface grid by inflating 
the previous surface at points that have continued to 
advance. 

10) Repeat steps 1 through 9 until no new points are 
accepted. 
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Fig. 2. Advancing-norma! point placement along 
normal line from boundary surface. 

Normal Spacing 

The normal spacing is determined using accelerated 
geometric growth. The initial normal spacing can be 
specified globally or at each boundary point. Standard 
geometric growth is used with an accelerated growth 
factor. The normal spacing is determined from 

As„+j =a„As„ (1) 

0^+/ = min (P a„ , (2) 

where As„ is the normal spacing for layer n, is the 
growth factor for layer n, is the maximum 

allowable growth factor, and p is the growth 
acceleration factor. 

Boundary Normal Vectors 

A boundary normal vector is required for the 
advancing-normal procedure. This normal is 
determined on the inflated surface during each pass. On 
the first pass the surface is the same as the original 
boundary surface. As only the boundary face normals 
are unique (since planar faces are used), some form or 
averaging or optimization procedure must be used to 
obtain the normal vector at nodes. Weighted averages 
can produce variations in normals due purely to 
topology or local face area differences. In the present 
work, a least-squares optimization procedure is used to 


eliminate those variations. An error function is defined 
as 

ej = I- bi • nj (3) 

where ej is the error function for face j, nj is the face unit 
normal vector for face y, bi is the node unit normal vector 
for node /. Node and face normals are shown in Fig. 3. 
The error function is also directly related to the volume 
of the element that will be produced from a given face. 
Minimizing the error function also maximizes the 
element volume. Least-squares optimization can be 
used to find bi such that S is minimized. The 
resulting equations are 

l[(bfnj)n/] = 'Ln/ (4) 

I [(bi • nj) njy] = Jjtjy (5) 

X f(bi • nj) n/J = I/i/ (6) 

where X denotes the sum over all faces surrounding 
node / and njy, and nj^ are the x, y, and z components 
of the unit normal vector for face j. 



Fig. 3. Node normals determined from surrounding 
face normals. 


Element Connectivity 

The element connectivity for the points created using 
advancing-normal point placement is determined 
directly from the order in which they were created. The 
initial point ordering is re-ordered so that optimal 
element quality will be produced (concave first and 
convex last). Tetrahedral elements are created (only 
temporarily if mixed elements are desired) for a given 
new point by inflating all surrounding boundary faces as 
shown in Fig. 4. Pentahedral elements can also be 
created using points generated on subsequent layers. 
Five and six node pentahedra are formed by combining 
tetrahedra as shown in Fig. 5. With mixed element 
types, five-node pentahedra and tetrahedra are created 
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only on the outer layer of the anisotropic region as 
shown in Fig. 6. These elements are required so that the 
anisotropic region is bounded only by triangular faces. 
All of the elements in the combined grid have strict 
node, edge, and face matching to each other and to 
neighboring tetrahedral elements. 



Fig. 4. Creation of tetrahedral elements using 
advancing-normal point placement. 



Fig. 5. Pentahedral elements formed by combining 
tetrahedral elements. 



Fig. 6. Use of pyramid and tetrahedral elements to 
transition between prismatic and tetrahedral 
regions. 


Isotropic Grid Generation 

The AFLR grid generation procedure used in the 
present work is a combination of automatic point 
creation, advancing type ideal point placement, and 
connectivity optimization schemes. A valid grid is 
maintained throughout the grid generation process. This 
provides a framework for implementing efficient local 
search operations using a simple data structure. It also 
provides a means for smoothly distributing the desired 
point spacing in the field using a point distribution 
function. This function is propagated through the field 
by interpolation from the boundary point spacing or by 
specified growth normal to the boundaries. Points arc 
generated using advancing-front type point placement. 
The connectivity for new points is initially obtained by 
direct subdivision of the elements that contain them. 
Connectivity is then optimized by local-reconnection 
with a combined Delaunay and min-max type 
(minimize the maximum angle) type criterion. The 
overall procedure is applied repetitively until a 
complete field grid is obtained. Complete details are 
presented in Refs. 7 and 8. 

A pplication Exampl es 

Selected application examples are presented here to 
demonstrate the capabilities of the present procedure for 
generation of three-dimensional unstructured grids of 
mixed element types that are suitable for 
Reynolds-Averaged Navier-Stokes simulations. All 
geometry preparation and surface grid generation work 
was done using SolidMesh^ with AFLR surface grid 
generation®. 

Grid quality distributions and statistics are presented for 
all examples in Fig. 7. Element angle is used as the grid 
quality measure. The complete set of grid quality data 
consists of the six, eight, and nine dihedral angles for all 
tetrahedra, five-node pentahedra, and six-node 
pentahedra respectively. In Fig. 7 the distribution plot is 
in 5 deg. increments. As shown, the distribution has 
peaks at 60 and 90 deg. from the six-node pentahedra 
and a peak near 70 deg. from the tetrahedra. The results 
for the examples presented are representative of those 
obtained for a variety of configurations. Typically, for 
isotropic elements in the grid, the maximum element 
angle is 160 deg. or less, the standard deviation is 17 
deg. or less, and 99.5% or more of the element angles are 
between 30 and 120 deg. And, for anisotropic elements 
in the grid, the maximum element angle is 170 deg. or 
less and 99.5% or more of the element angles are 
between 30 and 135 deg. The minimum angle is usually 
dictated by the geometry. Convex or concave edges with 
an included angle less than 20 deg., such as a sharp 
trailing edge or the interior of a wedge, can produce 
larger maximum angles in the anisotropic region. The 
maximum anisotropic element angle can be controlled 
by specifying the maximum allowable angle (which 
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eliminates generation of such elements) or by use of 
multiple normals at convex edgcs^. 



ELEMENT ANGLE 

Fig, 7, Volume element angle distributions. 


CPU time required on a SUN Ultra 60 workstation is 
presented in Table 1 for each example. Computer 
routines for the three-dimensional grid generator are 
written in C with dynamic memory that is automatically 
reallocated based upon actual requirements. All 
floating-point calculations arc performed using 64 bit 
precision with 8 byte data. The CPU times reported arc 
for one processor and include all I/O and generation of 
grid quality data. A boundary surface grid file is the 
input. The output includes a grid coordinate and 
connectivity file and a quality data file. Memory 
required is about 300 bytes per node generated. 
Requirements for memory and CPU time vary with the 
percentage of anisotropic elements, as those 
requirements for anisotropic generation are 
considerably less than those for isotropic generation. 

User input required to generate a complete grid is 
minimal and includes specifying the point spacing at 
selected control points on the boundary curves for 
surface grid generation. Selection of options such as 
which boundaries to generate anisotropic elements from 
and initial normal spacing are the only required user 
input for volume grid generation. There are no user 
adjustable parameters that need to be changed from case 
to case. In all cases presented here, the initial normal 
spacing was set suitable for high Reynolds number 
viscous CFD analysis. Initial normal spacing was 
determined such that the first node adjacent to a viscous 
surface would have a y'*' value near 1. 


Table 1. Number of nodes and elements generated and CPU time required for example cases. 


Case 

Boundary 

Faces 

Nodes 

Five-node 

Pentahedra 

Six-node 

Pentahedra 

Tetrahedra 

CPU Time 
(min) 

Destroyer 

Hull 

85,656 

1,269,048 

17,140 

2,234,663 

702,609 

9.6 

EETWing 

Body 

272,920 

2,290,661 

33,001 

3,744,105 

2,004,663 

41.9 

Space Shuttle 
Orbiter 

152,810 

1,102,869 

10,269 

1,709,948 

1,181,232 

16.6 

Titan IV-^A 
Exterior 

337,596 

2,980,493 

25370 

5,113,326 

1,962,705 

59.4 

Titan FV-A 
Interior 

88,502 

571,399 

31,548 

861,492 

642,859 

5.9 
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Destroyer Hull 

The first case presented is for the Navy model 5415 
destroyer hull. Multiple views of the surface grid on the 
water-line, hull, and propel lor surfaces are shown in 
Fig. 8. A field cut near the propellers is shown in Fig. 9. 
Element size varies smoothly in the field and there is a 
smooth transition between the anisotropic and isotropic 
regions. Grid quality distributions are shown in Fig. 7. 
Number of boundary faces, nodes, elements and CPU 
time are presented in Table I . 



b) stern and propellers 
Fig. 8. Destroyer hull surface grid. 



Fig. 9. Field cut near propellers for destroyer hull 
grid. 


Energy Efficient Tran sport (EET) 

A grid suitable for high Reynolds number viscous CFD 
analysis was generated for a high-lift wing/body 
configuration of the Energy Efficient Transport (EET). 
The surface grid on the underside of the wing is shown 
in Fig. 10. Field cuts near the wing/body, wing/slat, and 
wing/flap regions are shown in Fig. 11. Element size 
varies smoothly in the field and there is a smooth 
transition between the anisotropic and isotropic regions. 
In narrow regions between components, the number of 
anisotropic layers is reduced to produce high-quality 
elements between them. Also, the symmetry plane grid 
has been re-generated to match the interior anisotropic 
elements exactly. Grid quality distributions are shown 
in Fig. 7. Number of boundary faces, nodes, elements 
and CPU time are presented in Table 1. An all 
tetrahedral element version of this grid was used by 
Sheng, et al*^ for incompressible flow simulation with 
an implicit multi-block flow solver. 



Fig. 10. EET wing surface grid. 
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Fig. 11. Field cuts for EET grid. 



b) slat region 


NASA Space Shuttle Qrbiter 

A grid suitable for high Reynolds number viscous CFD 
analysis was generated for the NASA Space Shuttle 
Orbiter. The surface grid is shown in Fig. 12. Field cuts 
near the rocket motor and inboard flap regions are 
shown in Fig. 13, Element size varies smoothly in the 
field and there is a smooth transition between the 
anisotropic and isotropic regions. In narrow regions 
between components, the number of anisotropic layers 
is reduced to produce high-quality elements between 
them. Also, the symmetry plane grid has been 
re-generated to match the interior anisotropic elements 
exactly. Grid quality distributions are shown in Fig. 7. 
Number of boundary faces, nodes, elements and CPU 
time are presented in Table 1. 
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Fig. 12. NASA space shuttle orbiter surface grid. 



Fig. 13. Field cuts for NASA space shuttle orbiter 


a) rocket motor region grid. 
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Titan IV-A Launch Vehicle 

A grid suitable for high Reynolds number viscous CFD 
analysis was generated for a Titan IV-A launch vehicle 
wind tunnel test model configuration. This 
configuration includes two strap-on solid rocket motors 
(SRM), thrust vector control (TVC), stage separation 
motor (SSM), interstage cavity, and wind tunnel test 
model sting. The overall configuration and surface grid 
near the SRM, TVC and SSM regions are shown in Fig. 
14. Field cuts near the SRM, TVC, and SSM regions are 
shown in Fig. 15. A field cut within the interstage cavity 


is shown in Fig. 16. Element size varies smoothly in the 
field and there is a smooth transition between the 
anisotropic and isotropic regions. In narrow regions 
between components, the number of anisotropic layers 
is reduced to produce high-quality elements between 
them. Grid quality distributions are shown in Fig. 7. 
Number of boundary faces, nodes, elements and CPU 
lime are presented in Table 1. This configuration, and 
the others presented here, are representative of the high 
level of geometric complexity that can be handled 
routinely using the present approach. 



a) complete vehicle with wind tunnel sting 



c) TVC region Fig. 14. Titan IV-A launch vehicle surface and grid. 
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Summary 

A procedure has been presented for efficient generation 
of high-quality unstructured grids of mixed element 
types suitable for CFD simulation of high Reynolds 
number viscous flow fields. Layers of anisotropic 
elements are generated by advancing along prescribed 
normals from solid boundaries. The points arc 
generated such that either pentahedral or tetrahedral 
elements with an implied connectivity can be be directly 
recovered. As points are generated they are temporarily 
attached to a volume triangulation of the boundary 
points. This triangulation allows efficient local search 
algorithms to be used when checking merging layers. 
The existing AFLR procedure is used to generate 
isotropic elements outside of the anisotropic region. 
Results were presented for a variety of applications. The 
results demonstrate that high-quality anisotropic 
unstructured grids can be efficiently and consistently 
generated for complex configurations. 
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